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ABSTRACT 

An  abbreviated  review  of  rate  coefficients  relevant  to  HF  laser  kinetics  modeling  is  presented.  The  literature  has  been 
surveyed  from  the  last  published  review  in  1983  to  the  present.  Updated  HF  Einstein  emission  coefficients  are  tabulated. 
This  brief  summary  of  a more  detailed  review  addresses  rate  coefficients  relevant  to  HF  generation,  reactive  quenching, 
self-relaxation,  and  vibrational  relaxation  by  a selection  of  atoms  and  molecules.  In  addition,  a review  of  recent 
experiments  and  theoretical  calculations  relevant  to  the  role  of  rotational  non-equilibrium  in  HF  lasers  is  presented.  A list  of 
recommended  temperature  dependent  expressions  for  critical  reaction  rate  coefficients  is  given. 

1.  INTRODUCTION 

Since  its  invention  in  the  mid- 1960's,  the  HF  laser  system  has  been  extensively  studied  and  developed  to  the  point 
where  megawatt-class  devices  can  be  built.  In  fact,  most  of  the  research  in  the  recent  past  has  focused  on  large-scale  laser 
technology  demonstrations.  Despite  the  enormous  effort  expended  to  accomplish  this,  a complete  understanding  of  all 
facets  of  HF  laser  performance  is  still  evolving  and  is  not  complete.  For  example,  research  continues  into  the  role  of 
reagent  mixing  and  heat  transfer  between  the  fluids  and  the  construction  material  of  the  device.  Combustor  instabilities  and 
other  complex,  transient,  fluid  dynamical  features  also  impede  our  understanding  of  the  laser’s  performance. 

The  only  way  to  achieve  insight  into  the  details  of  the  HF  laser  is  to  employ  computational  fluid  dynamical  (CFD) 
codes  that  can  integrate  the  complex  fluid  properties  with  the  myriad  chemical  reactions  that  occur  in  the  laser  cavity. 
Unfortunately  (although  perhaps  not  surprisingly  considering  the  complexity  of  the  problem),  CFD  codes  have  had  limited 
success  at  accurately  modeling  real  HF  laser  systems.  As  a result,  both  the  laser  performance  data  and  the  reaction  rate 
constants  used  to  baseline  the  models  have  come  under  increased  scrutiny  in  recent  years.  This  scrutiny  has  uncovered 
serious  questions  about  the  kinetics  package  that  have  yet  to  be  answered  conclusively.  These  questions  include  the 
importance  of  rotational  nonequilibrium,  the  magnitude  of  various  quenching  processes,  the  role  of  three  body  and 
heterogeneous  fluorine  atom  recombination,  and  other  fundamental  properties  such  as  Einstein  coefficients. 

The  main  topics  of  this  report  (in  order  of  their  presentation)  are  Einstein  coefficients  and  relevant  kinetic 
measurements.  It  is  not  within  the  scope  of  this  document  to  discuss  fluid  dynamics  issues,  such  as  recently  developed  3 
dimensional  computational  fluid  dynamics  (CFD)  codes  or  new  algorithms  to  model  mixing  or  optical  resonators. 

2.  EXPERIMENTS  AND  CALCULATIONS  RELEVANT  TO  HF  LASER  MODELING 
2.1  Einstein  Coefficients 

The  Einstein  coefficients  used  by  most  HF  CFD  codes  are  based  on  the  values  found  in  the  Handbook  of  Chemical 
Lasersl(which  are  in  turn,  based  on  the  empirical  calculations  of  Herbelin  and  Emanuel2),  and  have  not  been  updated  in 
over  25  years.  Table  1 gives  a representative  sample  of  the  Handbook's  HF  vibration-rotational  Einstein  emission 
coefficients  as  well  as  the  more  recent  (and  preferred)  results  of  Setser  and  co-workers3,  see  below.  While  the  agreement  is 
generally  good  for  the  first  3 vibrational  levels,  large  differences  are  apparent  as  the  vibrational  quantum  number  increases. 
In  1991,  Zemke4  and  co-workers  published  a potential  surface  based  on  the  spectroscopically  determined  potential  of 
Coxon  and  Hajigeorgiou5,  adjusted  to  reproduce  the  proper  long-range  behavior  by  including  both  dispersion  and  exchange 
effects6.  In  the  same  publication,  Zemke  and  co-workers  provided  an  ab  initio  dipole  moment  function  that  spanned  the 
same  range  of  internuclear  distances  as  the  complete  potential  energy  curve.  The  resulting  Einstein  coefficients  should  be 
the  most  reliable  theoretical  values.  Shortly  after  the  publication  of  Zemke's  results,  Setser  and  co-workers  produced  an 
extensive  set  of  vibration-rotational  Einstein  A coefficients  for  HF/DF  and  HC1/DC13.  Their  calculations  used  an  RKR 
potential  and  the  ab  initio  dipole  moment  function  of  Ogilvie7.  Their  results  are  in  excellent  agreement  with  Zemke. 
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Considering  the  importance  of  the  Einstein  A coefficients  in  calculating  the  stimulated  emission  cross  section  and 
the  gain,  the  most  accurate  values  available  should  be  employed.  Unfortunately,  Zemke  and  co-workers  calculations  were 
only  for  a limited  number  of  rotational  quantum  numbers.  We  recommend  that  the  HF  CFD  codes  be  updated  with  the 
results  of  Arunan,  Setser,  and  Ogilvie3. 

2.2  HF  Kinetics 

Most  modern  day  HF  CFD  codes  have  kinetics  packages  that  are  based  on  a 1976  Aerospace  Corporation  technical 
report  by  N.  Cohen  & J.  Bott8.  This  report  and  its  1977  supplement9  contained  a thorough  review  of  contemporary 
literature  results  up  to  1977  and  recommendations  for  rate  constants  related  to  the  HF  laser.  An  update  was  published  5 
years  later  in  198210,  and  a few  relevant  reactions  were  reviewed  again  in  1983  by  Cohen  & Westberg11.  Beyond  these 
reports,  there  have  only  been  a handful  of  critical  evaluations  of  kinetic  data  relevant  to  the  HF  laser.  A 1982  review  article 
by  Leone  summarizes  hydrogen  halide  vibrational  energy  transfer  and  contains  rate  coefficients  relevant  to  the  HF  laser 
system12.  In  1983,  George  Hart  of  the  Naval  Research  Laboratory  reviewed  the  pulsed  DF  chemical  laser  codes  and  the 
corresponding  kinetic  database13.  Although  his  report  was  specifically  for  DF,  it  contains  a wide  variety  of  relevant  and 
helpful  evaluations  for  the  HF  laser.  The  following  paragraphs  summarize  the  most  recent  and  reliable  calculations  and 
experiments  relevant  to  HF  generation  and  quenching. 

2.2.7  HF  Generation  - H + F2  and  F + H2 

The  generation  of  HF(v)  in  the  HF  laser  can  proceed  via  one  of  two  reactions,  H + F2  or  F + H2,  which  have 
significantly  different  product  vibrational  distributions.  The  reaction  of  atomic  hydrogen  with  molecular  fluorine  (often 
referred  to  as  the  "hot"  HF  generation  reaction)  produces  highly  vibrationally  excited  HF,  while  F + H2  (the  "cold”  HF 
generation  reaction)  produces  only  moderate  vibrational  excitation,  see  below.  The  vibrational  distribution  for  H + F2 
peaks  at  v = 6 and  extends  up  to  v =9.  Table  3 summarizes  the  vibrational  distributions  recommended  by  the  Cohen  and 
Bott8-11  reviews  as  well  as  the  measured  distributions  from  a variety  of  experiments.  Most  experimental  measurements, 
particularly  those  of  Polanyi14,  Jonathan15,  and  Tardy16  analyzed  their  data  using  Einstein  coefficients  that  have  since  been 
shown  to  be  inaccurate2, 17, 18.  Hence,  corrected  distributions  using  the  recommended  set  of  A coefficients3  are  shown  in 
parentheses  in  Table  3. 

In  general,  the  available  experimental  results  are  in  reasonable  agreement  for  the  HF(v)  distribution.  The  only 
uncertainty  concerns  the  nascent  population  of  v > 8.  On  the  low  end,  the  fast  flow  reactor  studies  of  Setser19  and 
Kaufman20  found  no  Pg  - Pio,  while  on  the  high  end  the  pressure-pulse  chemiluminescence  mapping  experiments  of  Tardy16 
found  substantial  populations  for  v = 8 - 9.  The  presence  of  at  least  some  P8  - Pi0  is  supported  by  the  infrared 
chemiluminescence  studies  of  Polanyi  and  Jonathan,  who  reported  minor  P8  - P10.  The  nascent  vibrational  populations  from 
recent  theoretical  calculations21  are  in  satisfactory  agreement  with  experiment  but  have  slightly  narrower  distributions  with 
small  but  nonzero  population  of  v = 8 - 10. 

Surprisingly,  the  Cohen  and  Bott  reviews  recommend  no  initial  population  of  v = 0 - 2,  even  though  all  of  the 
experimental  measurements  (most  of  which  were  available  at  the  time)  indicate  small,  but  nonzero  V\  and  P2.  Clearly, 
some  initial  population  of  v = 1 - 2 is  indicated  by  the  experimental  evidence,  and  in  light  of  this,  we  recommend  the 
distribution  given  in  the  final  column  of  Table  2.  This  distribution  attempts  to  encompass  the  general  observation  that  v = 8 . 
- 9 is  present  but  at  lower  populations  than  suggested  by  Tardy.  The  recommended  values  for  v = i - 6 are  simply  the 
average  and  one  standard  deviation  from  the  5 experimental  measurements.  Recommendations  for  v = 0,  v = 7 - 10  are 
estimates  based  on  the  experimental  values  and  have  significantly  larger  error  bars  (±  50  % or  more).  Our  distribution  is 
similar  in  shape  to  the  DF(v)  distribution  generated  by  the  D + F2  reaction16, 22,  which  falls  off  rapidly  beyond  the  peak  at  v 
= 9-10. 

The  total  rate  constant  for  H + F2  has  not  been  firmly  established.  The  1982  Cohen  & Bott  review10  gives  k(T)  = 5.0  x 
10"i 5 T1,5exp(-845/T)  cm3  molecule'1  s'1  and  the  1983  Cohen  & Westberg  recommendation11  is  essentially  the  same,  4.8  x 
Iq-15  ji.4  eXp(_667/T)  cm3  molecules'1  s'\  In  both  cases,  the  recommended  value  was  based  on  the  experiments  of  Homann 
and  co-workers23  and  unpublished  transition  state  theory  calculations  of  Westberg  and  Cohen24.  The  1981  Baulch  kinetic25 
database  recommends  k = 1.46  x 10'10  exp(-I210fT)  cm3  molecule'1  s'1  for  T = 290  - 570  K and  points  out  that  Homann's 
result  is  significantly  smaller  (approximately  a factor  of  2 at  300  K)  than  previous  results  by  Rabideau26,  Vasil’ev27,  and 
Goldberg  . A new  measurement  of  the  total  H atom  removal  rate  constant  and  the  nascent  HF  distribution  was  performed 
recently  by  Heaven  and  co-workers29.  They  report  k = 2.4  ± 0.4  (2a)  x 10'12  cm3  molecule'1  s'1. 

The  F + H2  reaction  is  a prototypical  system  for  fundamental  reaction  dynamics,  and  as  such,  has  been  a favorite 
subject  for  both  theoretical  and  experimental  state- to- state  reactive  scattering  studies.  The  reaction  is  particularly  amenable 
to  molecular  beam  studies  and  vibrationally  state  resolved  differential  cross  sections  have  been  measured30'36. 
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Corresponding  high  level  ab  initio  calculations  and  simulations37, 38  have  achieved  very  good  agreement  with  experiment. 
The  vibrational  distribution  remains  unchanged  from  the  1982  Cohen  & Bott  report10, 0.00  : 0. 15  : 0.55  : 0.30  for  v = 0 - 3. 

A recent  review  by  Persky  & Komweitz39  has  refined  the  overall  rate  constant  for  the  F + H2  reaction.  Following  a 
detailed  examination  of  relevant  publications  they  recommend  k(T)  = 1.1  ± 0.1  x 10*10  exp(-(450  ± 50)/T)  cm3  molecule1  s' 
1 over  the  190  - 376  K temperature  range,  and  2.43  ± 0.15  x 10*n  cm3  molecule'1  s*1  at  298  K.  This  compares  reasonably 
well  with  literature  reviews  published  in  198311  (k(T)  = 4.5  x 10'12  T0  5 exp(-319/T)  and  k298  = 2.66  x 10'11  cm3  molecule'1  s' 
*),  199240,  and  199741  (k(T)  = 1.4  x 1010  exp(-(500  ± 200)/T)  and  k298  = 2.6  ± 0.6  x 1011  cm3  molecule'1  s'1),  as  well  as  the 
Cohen  & Bott  reviews  of  19779  (k(T)  = 2.7  x 10'10  exp(-805/T)  and  k298  = 1.8  x Iff11  cm3  molecule1  s'1)  and  198210  (k(T)  = 
4.32  x 10'12  T05  exp(-3G7/T)  and  k298  = 2.7  x 10'11  cm  molecule1  s'1).  The  most  recent  review  by  the  IUPAC 
Subcommittee  on  Gas  Kinetic  C>ata  Evaluation  for  Atmospheric  Chemistry42adopted  the  Persky  recommendation39. 
Unfortunately,  the  limited  temperature  range  of  the  Persky  expression  is  problematic  for  HF  laser  modeling,  since  the  laser 
typically  operates  at  substantially  higher  temperatures.  To  date,  there  has  been  only  one  experiment  that  has  measured  k(F 
+ H2)  above  376  K.  Heidner  and  co-workers43  monitored  the  time-resolved  infrared  emission  of  product  HF  following 
multi-photon  dissociation  of  SF6  in  the  presence  of  H2  over  the  295  - 765  K temperature  range.  The  resulting  Arrhenius 
expression  for  k(F  + H2)  is  2.2  ± 0.4  x 10'10  exp(-(595  ± 50)/T)  cm3  molecule*1  s'1,  just  7 % smaller  than  Persky's  at  room 
temperature  but  40%  larger  if  Persky's  expression  is  extrapolated  to  765  K.  Persky  & Komweitz39  considered  Heidner1  s 
results  "problematic  with  regard  to  the  calculated  kinetic  isotope  effect."  Indeed,  the  Heidner  experiment  gave  a 
temperature  independent  kinetic  isotope  effect  while  the  accepted  value  is 


kF+H  2 

- = 1.04  ± 0.02  exp((I86  ± 5)/T). 

kF+D2 


[1] 


In  lieu  of  more  data  for  T > 376,  we  recommend  the  conclusions  of  Persky  and  Komweitz39  for  the  190  - 376  K 
temperature  range  and  the  expression  of  Heidner  and  co-workers43  for  T > 376  K. 

Some  CFD  codes44  include  F atom  reactions  with  vibrationally  excited  H2  even  though  this  process  was  not  included  in 
the  original  Cohen  and  Bott  compilations.  There  have  been  no  specific  experimental  measurements  to  support  or  refute  this 
assumption  and  we  do  not  recommend  inclusion  of  reactive  processes  that  involve  vibrationally  hot  H2.  In  any  case,  it  is 
unlikely  that  inclusion  of  these  reactions  will  have  any  effect  on  the  overall  performance  of  the  laser  because  [H2(v  > 0)] 
should  be  extremely  small. 

2.2.2  Reactive  Quenching 

Vibrationally  excited  HF  can  be  removed  by  hydrogen  atoms  by  V-R,T  inelastic  collisions  or  by  chemical  reaction 
to  give  molecular  hydrogen  and  an  F atom: 

H + HF(v)  ->  H2  + F.  [2] 

In  principle,  microreversibility  enables  one  to  calculate  the  rate  constant  for  [2]  from  the  extensive  data  available  for 
the  well-studied  F + H2  reaction.  In  fact,  numerous  theoretical  studies  have  attempted  to  do  this  using  the  F + H2  potential 
energy  surface45-49.  According  to  these  calculations  the  barrier  to  F atom  transfer  is  large,  - 33  kcal  mol1,  and  reaction  [2] 
should  be  slow  for  v < 350'52.  This  is  consistent  with  the  experimental  results  of  Heidner  and  co-workers53'55  who  measured 
HF(v)  deactivation  by  H atoms  directly  using  HF  laser  induced  fluorescence,  and  the  flow  tube  measurements  of  Kwok  & 
Wilkins56.  These  experiments  report  a large  change  in  the  HF(v)  removal  rate  constant  for  v = 1-2  vs.  3.  This  change  is 
generally  attributed  to  the  opening  of  the  reactive  channel  for  v > 3.  However,  according  to  Heidner54, 55  only  a fraction  of 
the  total  H + HF(3)  encounters  that  result  in  removal  of  HF(3)  proceed  via  chemical  reaction  and  the  upper  limit  for 
reactive  quenching,  k(H  + HF(3)  ->  H2  + F),  is  5.0  x 10'11  cm3  molecule'1  s'1. 

The  1981  Baulch  kinetics  database25  makes  no  recommendation  for  k2  because  the  experimental  evidence  available  at 
the  time  was  inconsistent  with  the  data  for  the  well  established  forward  reaction,  F + H2.  While  there  have  been  no  new 
experiments  (for  thermal  collisions)  since  the  work  of  Heidner  and  Bott53'55,  the  available  theoretical  calculations  support 
their  slower  reaction  rate  constants8, 57. 


2.2.3  HF  Self-Relaxation 

One  of  the  most  active  areas  of  HF  kinetics  research  in  the  past  15  years  has  been  in  the  study  of  HF  self-relaxation  and 
vibrational  energy  transfer.  The  importance  of  these  processes  is  acute  because  this  is  the  dominant  relaxation  pathway  in 
the  HF  laser.  Unfortunately,  a consensus  regarding  the  magnitude  of  the  rate  constants  had  not  been  reached  prior  to  1982. 
For  example,  the  1977  Cohen  and  Bott9  compilation  contains  moderate  to  large  rate  constants  for  single-  and  multi- 
quantum V-R,T  deactivation  of  HF  by  ground  state  HF 

HF(v)  + HF  HF(v’)  + HF,  [3] 
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while  the  1982  Cohen  and  Bott  package10  includes  only  single  quantum  deactivations.  Significant  differences  also  exist  for 
the  HF  V-V  energy  transfer  reactions  such  as 

HF(v)  + HF(0)  4 HF(v-l)  + HF(1).  [4] 

Implicitly  included  in  reaction  [3]  are  the  V-R  redistribution  processes  that  produce  highly  rotationally  excited  HF: 

HF(v,  J)  + HF  HF(v-l,  J*  > 10)  + HF.  [5] 

As  will  be  discussed  in  detail  below,  this  quenching  process  is  of  special  interest  because  it  has  the  potential  to  contribute 
significantly  to  rotational  nonequilibrium. 

Shortly  after  the  publication  of  the  1982  review10,  the  Crim  group  at  the  University  of  Wisconsin58'62  and  the  Kaufman 
group  at  the  University  of  Pittsburgh20,63-66  undertook  a major  effort  to  characterize  the  total  self-relaxation  rate  constants 
and  the  mechanism  for  HF  self^relaxation.  Crim's  group  used  a double  resonance  type  of  experiment  where  the 
vibrationally  excited  HF  molecules  were  prepared  in  discrete  ro- vibrational  states  by  a pulsed  laser.  The  time  resolved 
fluorescence  and/or  Av  = 1 absorption  signals  were  analyzed  to  determine  total  quenching  rates  and  relaxation  mechanisms. 
Kaufman,  on  the  other  hand,  prepared  vibrationally  excited  HF  in  a flow  reactor  where  dilute  flows  of  H or  F atoms 
(generated  by  a microwave  discharge)  were  reacted  with  a variety  of  F or  H atom  donors.  The  IR  emission  was  collected 
with  an  InSb  detector  and  circularly  variable  filter.  A modified  Stem- Volmer  analysis  was  applied  to  the  quenching  data. 
Table  3 compares  the  experimentally  determined  rate  constants  for  HF  self-relaxation  with  a variety  of  other  experiments67' 
78,  relevant  calculations79'83,  and  the  standard  kinetics  packages8'11.  The  agreement  for  v = 1-7  is,  in  general,  excellent  and 
k3  is  well  established.  The  experimental  relaxation  rates  scale  as  v2'9  and  are  independent  of  the  initial  rotational  quantum 
number. 

In  addition  to  total  quenching  rate  constants,  the  Crim  and  Kaufman  laboratories  also  determined  the  relaxation 
mechanism.  Kaufman's  group  argued  strongly  for  a V-T,R  mechanism  (rather  than  V-V  energy  transfer)  based  on 
Lambert-Salter  plots20, 63-66  and  the  magnitude  of  the  rate  constants.  In  particular,  they  pointed  out  that  if  the  predominant 
mechanism  were  V-V  energy  transfer,  (eg.  HF(7)  + HF(0)  ->  HF(6)  + HF(1))  then  the  rate  constant  for  the  exothermic 
reverse  process  would  be  100  times  greater  than  the  gas  kinetic  limit.  Crim's  double  resonance  experiments  were  able  to 
quantify  the  role  of  V-V  energy  transfer.  They  found  that  the  fraction  of  inelastic  HF(v)  + HF(0)  room  temperature 
encounters  thatproceed  via  V-T,R  relaxation,  is  1.0,  0.41  ± 0.10, 0.56  ± 0.05, 0.84  ± 0.05,  and  0.98  ± 0.19  for  v = 1 - 5, 
respectively59,  ’ 62.  For  v > 3,  vibrational  energy  transfer  to  the  ground  state  collision  partner  plays  a relatively  minor  role 
in  the  relaxation  process.  Both  Crim  and  Kaufman  agree  that  multi-quantum  relaxation  is  unimportant62, 63  even  though 
work  by  Pimentel  and  Thompson  (see  below)  suggested  the  possibility  of  multi-quantum  V-R  transfer  with  Av  as  large  as  5. 
Crim  probed  the  role  of  multi-quantum  deactivations  directly  and  found  that  0.98  ± 0.19  and  0.87  ± 0.21  of  the  relaxed 
HF(4)  and  HF(5)  molecules,  respectively,  appear  in  the  next  lower  vibrational  level62.  It  is  important  to  note  that  Crim  and 
co-workers'  results  are  based  on  the  assumption  that  V-T,R  processes  that  produce  metastable  high  rotational  states  (which 
would  not  be  detected  in  their  experiment)  can  be  neglected.  The  invariance  of  the  vibrational  relaxation  rate  constant  with 
initial  rotational  quantum  number  and  the  work  of  Leone  (see  below)  tend  to  validate  this  assumption. 

Finally,  Crim  and  co-workers  found  that  the  rate  constants  for  HF  self-relaxation  are  inversely  dependent  on 
temperature59, 61 . Crim  and  co-workers  interpret  their  temperature  dependent  data  in  terms  of  relaxation  probabilities.  The 
functional  form  of  the  fitting  function  suggests  that  long-range  forces  dominate  the  relaxation  process: 

Pv(T)  = AT~m  [6] 

where  Pv  = kjkc  (kc  is  the  gas  kinetic  rate  constant  for  a collision  diameter  of  0.25  nm)  and  A and  m are  fitting  parameters. 
Strangely,  the  A values  determined  by  our  fits  (where  m = 1.3  was  fixed)  vary  significantly  from  Crim's  analysis61: 

A(Crim)  = 22,  370,  880,  and  1850  for  v = 1,  3, 4,  and  5 respectively,  while  A(this  work)  = 12,  315,  764,  and  1610  for  the 
same  v levels.  Nonetheless,  considering  the  overall  agreement  in  the  literature  for  v = 1 - 7 and  the  accuracy  of  double 
resonance  technique,  we  recommend  Crim's  HF  self-quenching  rate  constants61.  The  temperature  dependence  of  v > 6 has 
not  been  measured.  If  the  temperature  dependence  found  for  v = 1 - 5 (i.e.  kv  = Pv*kc  = kc*A*T1’3)  is  applied,  then  A(6)  = 
3107  and  A(7)  = 4339  are  calculated  from  the  measured  room  temperature  values20. 

2.2.4  HF  Relaxation  and  V-V  Energy  Transfer  with  H2 

There  have  been  several  studies  of  HF(v)  relaxation  by  a variety  of  molecular  quenchers.  Table  4 summarizes  the 
results  for  Q = H2,  which  are  the  most  relevant  to  HF  laser  kinetics20, 7I* 75, 84’87.  With  the  exception  of  the  work  by  Poole 
and  Smith75,  the  agreement  for  the  total  quenching  rate  constants  is  good  for  v = 3 - 5.  The  mechanism  for  the  quenching  is 
generally  believed  to  be  V-T,R  for  v = 3 - 5 for  two  important  reasons.  First,  V-V  energy  transfer  from  HF(v)  to  H2  is 
endothermic  for  all  single  vibrational  quantum  changes  in  HF: 
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HF(v)  + H2(0)  + AE  > 200  cm'-»  HF(v-l)  + H2(l)  [7] 

In  fact,  due  to  the  anhamonicity  of  HF,  the  energy  gap  between  HF(Av  = -1)  and  H2(0-l)  increases  with  vibrational 
quantum  number  (AEv.v  = -198  cm’1  for  HF(1)  and  *1171  cm’1  for  HF(7)).  Secondly,  if  the  endothermic  V-V  process  were 
the  dominant  mechanism,  then  the  exothermic  reverse  process 

HF(v-l)  + H2(l)  -»  HF(v)  + H2(0)  [8] 

would  be  several  times  larger  than  the  gas  kinetic  limit.  Only  for  HF(1)  does  a V-V  process  seem  possible,  and  indeed, 
vibrational  energy  transfer  is  the  most  likely  mechanism  for  v = 1.  For  v > 1,  however,  the  V-TJR  process 

HF(v)  + H2“»HF(v-l)  +H2  [9] 

should  be  the  dominant  mechanism. 

The  1982  Cohen  and  Bott  review10  contains  temperature  and  vibration  dependent  expressions  for  reactions  [13]  and 
[15]  even  though  there  is  no  convincing  evidence  for  a temperature  dependent  quenching  process85. 

k1{refA\)=2Ax\QXQvQ35TQSe(A01'li£v)IRT  cm^morXs~X  [10] 

k9(ref.  11)  = v27(o.6xl012r_1  +1.0x1 04r228)cmW_15_I  [11] 

For  example,  in  1973  Cohen  & Bott  measured  the  temperature  dependence71  (T  = 295, 450  - 1000  K),  of  HF(1)  deactivation 
by  H2and  found  that  the  total  deactivation  (k7  +k9)  rate  was  independent  of  temperature  (see  Figure  4 of  ref 71).  A year  later 
Bott  re-measured  the  temperature  dependent  quenching  of  HF(1)  from  440  - 690  K and  the  data  showed  considerable 
scatter  and  only  a weak  temperature  dependence84.  Finally,  in  1980  Bott  and  Heidner  measured  HF(1)  and  HF(3) 
quenching  by  H2  at  295  and  200  K and  found  deactivation  rate  coefficients  that  were  constant  vs.  T within  their 
experimental  error85.  Clearly,  there  is  not  sufficient  evidence  to  support  a T dependent  quenching  rate  constant. 

Cohen  & Bott’s  expression9  for  k9  significantly  underestimates  the  measured  values  at  room  temperature,  see  Table  4. 

In  fact,  the  1982  package10  eliminates  the  energy  transfer  reactions  from  H2(2)  and  H2(3)  as  well  as  the  Av  > 1 exchanges 
contained  in  the  1977  Cohen  and  Bott  package  because  there  is  no  specific  experimental  justification  for  them.  We 
recommend  the  rate  constant  values  listed  in  Table  4 for  the  V-T,R  quenching  of  HF  by  H2  and  assign  a T°  temperature 
dependence.  The  HF(v)  + H2(v')  V-V  energy  transfer  reactions  (reaction  [8])  are  calculated  from  detailed  balance. 

23  Rotational  Non-equilibrium 

The  question  of  rotational  non-equilibrium  for  the  HF  laser  system  has  been  the  subject  of  considerable  controversy  for 
many  years88, 89.  The  presence  of  rotational  non-equilibrium  in  the  HF  laser  was  first  suggested  by  pulse  initiated  HF  laser 
experiments  by  Pimentel  and  co-workers90*93  which  generated  lasing  on  HF  rotational  transitions  with  J as  high  as  33  in  the 
v = 1 manifold  and  J = 29  in  the  v = 0 manifold.  The  observation  of  rotational  laser  emission  is  an  extremely  sensitive 
method  for  studying  rotational  occupancies  because  the  population  inversions  needed  to  produce  the  laser  emission  are  100 
times  lower  than  for  ro-vibrational  transitions93.  Pimentel's  analysis  of  the  transient  behavior  of  the  laser  emission 
suggested  that  collisional  V-T,R  energy  transfer  reactions  that  populate  the  high  rotational  states  were  responsible  for  the 
observed  positive  gain,  rather  than  direct  pumping  by  the  initiating  reaction.  A remarkably  similar  phenomenon  was 
observed  by  Robinson  and  co-workers  in  their  work  with  HF94,  OH95,  and  NH%.  Further  experimental  observations  of 
emission  from  high  J states  following  HF(v)  quenching  by  CO,  C02,  and  HCN97, 98  as  well  as  quasiclassical  trajectory 
calculations99, 100  also  support  the  assertion  that  high  rotational  states  are  produced  in  the  V-T,R  relaxation  process. 

The  evidence  is  clear  that  the  principal  HF(v)  relaxation  mechanism  in  the  HF  laser  environment  is  HF  self-relaxation, 
and  that  the  relaxation  proceeds  via  V-T,R  energy  transfer.  There  are  two  important  questions,  however,  that  remain: 

1)  What  are  the  specific  products  of  the  V-T,R  process?  Figure  I summarizes  the  possible  relaxation  / energy  transfer 
routes  for  HF(v  = 2).  The  possible  mechanisms  include  "true"  V-T,R  relaxation  (the  solid  arrow)  where  the  loss  of  a 
vibrational  quantum  results  in  some  small  amount  of  rotational  and  translational  energy  transfer  to  the  HF(v  = 0)  quencher 
or  near-resonant  V-R  redistribution(the  broken  arrows),  where  the  quenched  HF  molecule  relaxes  to  a lower  vibrational 
state  with  a high  rotational  quantum  number  and  very  little  energy  is  transferred  to  the  quencher. 

2)  Are  the  high-J  HF  molecules  produced  by  the  self  -relaxation  process  "metastable”?  In  general,  rotational  relaxation 
rate  constants  (k  ~ 10*10  - I0'9  cm3  molecule'1  s’1)  are  10  - 100  times  larger  than  vibrational  deactivation  rate  constants  (k  - 
10'12  - 10’10  cm3  molecule'1  s’1).  However,  because  the  separation  between  HF  rotational  levels  is  large,  it  is  possible  that 
kR.R/r  ~ kv-R.T  for  sufficiently  high  J levels.  If  so,  the  vibrational  relaxation  process  could  significantly  perturb  the 
equilibrium  rotational  distribution  and  considerable  errors  could  be  realized  when  attempting  to  model  real  HF  laser 
devices. 

2.3.1  Relevant  Experimental  Studies 

As  was  discussed  above,  the  self-relaxation  measurements  of  Crim62,  Kaufman63,  and  Moore77  all  concluded  that  multi- 
quantum deactivations  were  not  important.  In  particular,  Crim  determined  that  0.98  ± 0.19  and  0.87  ± 0.21  of  the  relaxed 
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HF(4)  and  HF(5)  molecules,  respectively,  appear  in  the  next  lower  vibrational  level62,  Kaufman  came  to  the  same 
conclusion  and  suggested  that  multiquantum  relaxation  processes  account  for  less  than  1%  of  the  total  measured  relaxation 
rate  constants63.  Thus,  any  high-J  states  that  are  produced  by  the  HF(v)  self-relaxation  process  will  almost  certainly  be 
found  in  the  J = 0 - 20  range  of  the  next  lower  vibrational  level. 

The  role  of  V-TJR.  relaxation  reactions  that  populate  high  rotational  states  of  the  next  lower  vibrational  level  was 
addressed  directly  by  the  work  of  Haugen,  Pence,  and  Leone101  who  measured  the  time  dependent  population  of  HF(v  = 0,  J 
= 10  - 14)  following  pulsed  generation  of  HF(v  = 1,  J = 6).  They  concluded  that  a substantial  fraction  of  the  relaxation  of  v 
= 1 occurs  through  the  high  lying  rotational  levels  of  v = 0 (-20  - 40%  of  the  total  v = 1 relaxation  rate).  The  total 
phenomenological  self-relaxation  rate  constant  for  HF(v  = 1)  (which  by  definition  for  v = 1 is  purely  V-T,R)  that  they 
measured  was  identical  to  that  determined  in  the  double  resonance  experiments  of  Crim  & co-workers58,  (k  = 1.46  ± 0.1  x 
10*12  cm3  molecule*1  s’1). 

There  has  been  a significant  effort  in  the  last  20  years  to  measure  and  predict  rotational  relaxation  rate  constants.  Most 
recently,  Muyskens,  Copeland,  and  Crim102  106  have  measured  rotational  relaxation  rate  constants  for  HF(v  = 2 - 4,  J = 0 - 4) 
with  a variety  of  colliders.  Their  results  generally  confirm  the  standard  view  that  rotational  relaxation  is  10  -100  times 
faster  than  vibrational  relaxation,  particularly  for  the  lower  rotational  quantum  numbers  (J  = 0 - 8).  In  addition  to  probing 
the  role  of  V-T,R  relaxation,  the  Leone  group101  has  also  measured  rotational  relaxation  rate  constants  for  v = 0,  J = 10  - 14. 
They  found  no  experimental  evidence  of  bottleneck  effects  and  concluded  that  R-R,T  rates  always  exceed  the  V-TJR  rate  by 
one  or  two  orders  of  magnitude,  even  for  J = 10  - 13.  While  their  initial  report101  recommended  R-R,T  rate  constants  for  v 
ss  0,  J = 10  - 14  which  range  from  -1.2  x 10*10  - 6 x 10*11  cm3  molecule*1  s*1,  subsequent  measurements  in  the  same 
laboratory  suggested  even  larger  values107, 108.  The  dominance  of  R-R,T  relaxation  over  V-TJR  and  V-V  energy  transfer 
extends  to  other  colliders  besides  HF.  Taatjes  and  Leone,  for  example,  measured  the  rotational  relaxation  rate  constants  for 
HF  with  a variety  of  collision  partners  (Ar,  He,  Ne,  Kr,  Xe,  H2,  and  D2)107  and  found  that  while  rotational  relaxation  by 
atomic  species  is  very  inefficient  relative  to  HF,  H2,  and  D2  the  rotational  relaxation  rate  constants  for  atomic  quenchers 
exceed  the  vibrational  deactivation  rate  constants  by  several  orders  of  magnitude:  k(R-R,T)  > 10‘12  and  k(V-R,T)  = 10*17  - 
10*18  cm3  molecules*1  s*1.  Leone  and  co-workers  extended  their  measurements  to  non-ambient  temperatures108  and  found 
that  the  HF  V-T,R  and  R-R,T  self  relaxation  reactions  have  a negative  temperature  dependence.  The  negative  temperature 
dependence  for  rotational  relaxation  of  HF(v  = 0,  J = 13)  is  dramatic,  T 15  . 

In  addition  to  the  direct  experimental  measurements,  there  have  also  been  attempts  to  extrapolate  the  low-J  results  to 
high-J  using  scaling  laws  and  approximations  such  as  the  exponential  energy  gap  law(EEG),  the  power  law  model  (PLM), 
and  the  energy  corrected  sudden  (ECS)  approximation.  Most  of  these  efforts  are  summarized  elsewhere88, 89,  and  while  the 
accuracy  of  the  models  for  predicting  accurate  R-R,T  rates  is  the  subject  of  some  controversy,  two  general  conclusions  may 
be  drawn  from  the  relevant  literature.  1)  The  PLM  and  ECS  models  give  the  most  reliable  results  when  compared  to  the 
existing  high  J and  low  J data.  The  EEG  model  consistently  underestimates  k(R-R,T),  in  some  cases  by  several  orders  of 
magnitude.  2)  The  rate  constants  for  rotational  relaxation,  k(R-R,T),  are  large,  > 10'11  cm3  moleculs*1  s*1. 

One  noteworthy  pair  of  reports109,  uo,  which  claim  to  use  a "more  reliable  form  of  the  power  scaling  law"  to  calculate 
rotational  energy  transfer  rate  constants  for  v = 1 - 2,  J = 0 - 20  give  k(R-RT)  values  on  the  order  of  10*11  - 10*10  cm3 
molecule*1  s"\  even  for  J = 20.  On  the  other  hand,  their  results  suggest  that  rotational  relaxation  rates  actually  increase  with 
vibrational  energy,  contrary  to  the  results  of  Crim  and  coworkers10  *106.  The  reliability  of  their  model  is,  as  the  authors 
themselves  admit,  "an  open  discussion." 

2.3,2  Relevant  Theoretical  Studies 

No  review  of  the  role  of  rotational  equilibrium  for  the  HF  laser  would  be  complete  without  some  discussion  of 
quasiclassical  trajectory  calculation  results,  most  notably  those  of  Wilkins  and  Kwok79, 80,  m*  n\  Thompson99, 100, 113, 114,  and 
Billing83,  ll5'U7.  Billing’s  calculations  found  no  evidence  of  high  rotational  state  population,  while  the  calculations  by 
Wilkins  and  Thompson  indicate  that  vibrational  -rotational  energy  transfer  is  a relatively  efficient  process  and  that  multi- 
quantum deactivations  occur  on  a fairly  regular  basis.  In  particular,  Thompson99, 100  calculated  state-to-state  collsion  cross 
sections  for  HF(v  = 4,  J =20)  relaxation  by  He  and  reported  3.3,  6.7,  10.7,  18.5,  and  38.24  a.u.2  for  Av  = 4,  3,  2,  1 and  0, 
respectively.  Calculations  of  this  sort  are  usually  very  sensitive  to  the  details  of  the  potential  surface  on  which  the 
trajectories  run,  and  unfortunately  the  requisite  state-to-state  cross  sections  required  to  evaluate  the  reliability  of  the 
theoretical  calculations  have  yet  to  be  measured.  In  general,  the  available  experimental  data  does  not  support  multi- 
quantum deactivations. 

In  summary,  the  majority  of  the  available  evidence  supports  single  vibrational  quantum  V-T,R  relaxation  which 
populates  the  high  rotational  states  of  the  next  lower  vibrational  state.  There  is  no  specific  experimental  evidence 
supporting  multiquantum  vibrational  V-T,R  relaxation.  There  is  no  doubt  that  near-resonant  V-T,R  relaxation  plays  an 
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important  role  in  the  HF  chemical  laser  system  and  successful  quantitative  modeling  depends  on  its  inclusion.  However, 
the  available  experimental  evidence  clearly  shows  that  k(V-TJR)  « k(R-R,T),  and  in  light  of  this,  it  is  doubtful  that  high  J 
states  can  act  as  "reservoirs”  for  near-resonant  lasing  levels.  It  seems  more  likely  that  the  V-T,R  process  simply  reduces  the 
gain  of  the  (1-0)  and  (2-1)  transitions  by  reducing  the  population  of  the  upper  state  while  simultaneously  increasing  the 
population  of  the  lower  state. 

3.  CONCLUSIONS 

Table  5 summarizes  the  recommendations  of  this  report.  Overall,  many  of  the  expressions  found  in  the  1982  review  by 
Cohen  & Bott10  remain  valid  today,  in  particular,  the  elimination  of  multi-quantum  deactivation  reactions  that  were  a key 
feature  of  the  1977  kinetics  package.  These  kinds  of  relaxation  processes  have  been  demonstrated  to  be  very  slow  and  can 
be  safely  neglected.  Other  areqf  of  agreement  include  the  total  HF  generation  rate  constants  and  the  relaxation  rate 
constants  for  collisions  with  molecular  and  atomic  quenchers.  A new  measurement  of  the  H atom  removal  rate  constant  for 
the  "hot"  reaction,  H + F2,  would  be  particularly  useful. 

The  major  changes  that  we  suggest  occur  in  the  Einstein  coefficients,  HF  self-relaxation,  and  the  nascent  distribution 
for  H + F2.  While  in  many  cases  these  changes  are  minor,  they  may  ultimately  have  significant  effects  to  CFD  calculation 
results  due  to  enormous  complexity  of  the  HF  laser  system. 

Clearly,  there  are  some  aspects  of  the  HF  kinetics  package  that  should  be  re-examined  experimentally.  For  example,  in 
the  case  of  HF(v)  + F,  H,  Ar,  and  He,  the  recommended  expressions  are  based  on  only  a handful  of  measurements  at  a 
narrow  range  of  temperatures.  While  the  role  of  multi-quantum  deactivations  is  very  small  according  to  the  available 
experimental  data,  some  believe44  that  the  v19  scaling  law  for  the  HF  self-relaxation  process  may  be  indicative  of  open 
multi-quantum  deactivation  relaxation  pathways,  particularly  for  high  v.  Direct  measurements  for  the  Treanor  pumping 
(reaction  [14])  rate  constants  are  also  needed,  particularly  for  HF(v  > 1)  + HF(v  > 1),  for  which  no  data  currently  exists. 
Clarification  of  these  issues  would  undoubtedly  significantly  enhance  our  understanding  of  the  HF  laser. 
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Table  2:  Experimentally 
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<0.04  (<0.03) 

<0. 10  (<0.08) 

0.00  (0.00) 

0.00 

0.00 

0.04  ± 0.04 

1 

0.00 

0.00 

0.09  (0.06) 

0.12(0.07) 

0.15(0.14) 

0.07 

0.06 

0.08  ± 0.03 

2 

0.00 

0.00 

0.11  (0.08) 

0.13(0.10) 

0.13(0.12) 

0.17 

0.12 

0.13  ±0.03 

3 

0.18 

0.21 

0.13(0.10) 

0.25  (0.20) 

0.27  (0.26) 

0.28 

0.17 

0.20  ±0.07 

4 

0.30 

0.39 

0.45  (0.36) 

0.35  (0.30) 

0.41  (0.40) 

0.59 

0.37 

0.40  ±0.11 

5 

0.80 

0.70 

0.89  (0.83) 

0.78  (0.70) 

0.72  (0.70) 

0.93 

0.76 

0.78  ±0.10 

6 

1.00 

1.00 

1.00(1.00) 

1.00(1.00) 

1.00(1.00) 

1.00 

1.00 

1.00 

7 

0.00 

0.45 

0.45  (0.43) 

0.40  (0.48) 

0.76  (0.80) 

0.52 

0.62 

0.50  ±0.25 

8 

0.00 

0.36 

0.20(0.19) 

0.26  (0.37) 

0.46  (0.49) 

0.00 

0.00 

0.30  ±0.15 
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0.00 

<0.04  (<0.01) 

0.16(0.12) 

0.41  (0.43) 

0.00 

0.00 

0.15  ±0.15 

10 

0.00 

0.00 

<0.04  (<0.01) 

0.00  (0.00) 

0.00  (0.00) 

0.00 

0.00 

0.01  ± 0.01 

a The  values  in  parentheses  for  Jonathan,  Polanyi  and  Tardy  are  corrected  for  the  Einstein  coefficients  of  Setser3. 
b The  corrected  values  shown  were  calculated  from  the  distributions  reported  by  Kaufman20  which  were  corrected  for 
the  Einstein  coefficients  of  Sileo  & Cool118. 

c Tardy16  originally  used  the  Einstein  coefficients  of  Meredith  and  Smith17. 
d See  text  for  details 
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Table  3:  HF  total  self  relaxation  rate  constants 


Reference 

k3oo(HF(v)  + HF)  products,  (10‘12  cm3  molecules'1  s'1) 

v = 1 

v = 2 

v = 3 

v = 4 

v = 5 

v = 6 

WB3MI 

Experiments  f 

Bott  & Cohen'1 

1.8  ±0.3 

Hinchen  & Hobbs69 

1.8  ±0.2 

Bina  & Jones70 

2.3  ±0.3 

5 ± 2 

Kwok  & Wilkins72 

1.6  ±0.6 

16  ±5 

26  ±9 

27  ±10 

(58)a 

(loiy 

Osgood,  et  al.76 

1.7 

25  ±7 

49  ± 15 

43  ± 18 

Airey  & Smith73 

16 

17 

>44 

>65 

Poole  & Smith74’ 75 

13 

19 

32 

46 

52 

-43 

77 

Douglas  & Moore 

28  ±4 

72  ±5 

Lampert  et  al.78 

32  ±6 

88  ±11 

v , 20, 63-66 

Kaufman 

1.8 

19 

31 

73 

140 

290 

450 

Sft 

Copeland,  et  al. 

1.46  ±0.1 

19.8  ± 1.0 

Jursich  & Crim62 

30.2  ± 3.0 

72.8  ± 2.7 

151  ±8 

Calculations  | 

Wilkins  & Kwok/y’ 80 

1.7 

22 

29 

33 

42 

51 

Coltrin  & Marcus81, 82 

0.2  ±0.1 

19  ±3 

28  ±4 

53  ±10 

69  ± 10 

156  ± 11 

455  ±49 

Billing  & Poulsen83 

0.81 

6.2  ±2.2 

10  ±4 

19  ±7 

27  ±10 

43  ± 15 

82  ±29 

Standard  Kinetics  Packages1*  1 

Cohen  & Bott  19779 

MEM 

6.62 

9.94 

4.97 

16.6 

23.2 

82.8 

Cohen  & Bott  198210 

MEM 

10.0 

28.8 

60.9 

108.8 

174.7 

260.9 

1 1 M.  A.  Kwok  and  N.  Cohen,  personal  communication  reported  in  M. 

j b only  single  quantum  deactivation  rate  constants  are  listed. 

Table  4:  Room  Temperature  Quenching  Rate  Constants  for  HF  + H2 


Reference 

k100(HF(v)  + H2  ->  HF(v-l)  + H2  (1012 

cm3  molecules'1  s'1) 

v = 1 

v = 2 

v = 3 

v = 4 

v = 5 

< 

II 

o\ 

v = 7 

Bott  & Cohen7*’ 84 
Poole  & Smith75 

0.52  ± 0.03 

0.21 

0.15 

0.21 

0.49 

0.99 

1.6 

Douglas  and  Moore86 

0.31  ± .06 

0.47  ±.12 

Bott  & Heidner85 
Kaufman20 

0.52  ±0.05 

0.35  ± .04 

1.7  ±0.5 

3.5  ± 1 

9.1  ±2.7 

87 

Jursich,  et  al. 

0.38  ± 0.25 

0.67  ±0.10 

1.64  ±0.19 

0.01 

0.07 

0.21 

0.45 

0.83 

1.35 

Recommended 

0.52  ±0.04 

0.20  ±0.04 

0.35  ± 0.04 

0.50  ±0.2 

1.6  ±0.3 

3.5  ±1 

9.1  ±2.7 
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Figure  Is  Detailed  V-T,R  relaxation  pathways.  The  distinction  between  "true"  V-T,R  energy  transfer  and  V-R 
redistribution  is  shown  by  the  solid  and  broken  lines,  respectively.  Numerous  combinations  of  HF(1,  J)  and 
HF(0,  J)  states  can  be  populated  by  HF  V-T,R  energy  transfer.  Because  of  the  presence  of  near-resonant  energy 
levels  in  v = 1 and  0,  V-R  redistribution  can  populate  high  rotational  states  of  v = 0 and  1.  For  example,  the  near 
resonant  V-R  redistribution  pathways  shown  in  the  figure  have  energy  defects  of  -46.9,  519.6,  and  273.7  cm’1, 
for  relaxation  to  (1,15),  (1,14)  and  (0,20),  respectively.  On  the  other  hand,  if  HF(2,  6)  is  relaxed  to  HF(1,  6)  by 
HF(0,  J),  up  to  8 quanta  of  rotational  energy  can  be  transferred  to  the  HF(v  = 0)  molecule. 


• J = 7 


t j r 

0 1 2 

vibrational  quantum  number 
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